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Abstract 

The concept of the "transient effect" is examined in respect of a "mean first 
passage time". It is demonstrated that the time the fissioning system stays inside 
the barrier is much larger than suggested by the transient time, and that no en- 
hancement of emission of neutrons over that given by Kramers' rate formula ought 
to be considered. 
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1 Introduction 

Fission at finite thermal excitation is characterized by the evaporation of light particles 
and 7's. Any description of such a process must rely on statistical concepts, both with 
respect to fission itself as well as with respect to particle emission. For decades it has 
been customary to describe experiments in terms of particleQ and fission widths, where 
the former, F^, is identified through the evaporation rate and the latter Fj is given by 
the Bohr- Wheeler formula Fj = Fbw for the fission rate. Often in the literature this is 
referred to as the "statistical model". It was only in the 80's that discrepancies of this 
procedure with experimental evidence was encountered: Sizably more neutrons were seen 
to accompany fission events than given by the ratio F„/Fbw (for a review see e.g. [0]). It 
was then that one recalled Kramers' old objection against the Bohr- Wheeler formula. 
Indeed, in this seminal paper he pointed to the deficiency of the picture of Bohr and 
Wheeler in that it discards the influence of couplings of the fission mode to the nucleonic 
degrees of freedom. Such couplings will in general reduce the flux across the barrier as the 
energy in the fission degree of freedom Q may be diminished and fall below the barrier. 

^To simplify the discussion we will not distinguish the nature of the "particles" and in this sense 
include 7's in this notation. 



In Kramers' picture this effect is realized through the presence of frictional and fluctu- 
ating forces (intimately connected to each other by the fluctuation dissipation theorem). 
For not too weak dissipation Kramers' rate formula writes 

Tk = ^ -P ) (\A^ - = ^ = ««K. (1) 

Here, T and E-^, stand for temperature and barrier height, Wa for the frequency of the 
motion in the (only) minimum at Q = <5a and ?7b = {'~f /{2Mw))h for the dissipation 
strength at the barrier (at Q = Qb) with 7 being the friction coefficient and M the inertia. 
For the sake of simplicity we will assume these coefficients not to vary along the fission 
path; otherwise the formula must be modified 0. For vanishing dissipation strength (|I]) 
reduces to the Bohr- Wheeler formula (simplified to the case that the equilibrium of the 
nucleons can be parameterized by a temperature). 

Commonly, formula (0) is derived in a time dependent picture solving the underlying 
Fokker-Planck equation for special initial conditions with respect to the time dependence 
of the distribution function^. The initial condition is intimately related to the condition of 
a compound reaction, in that the decay process is assumed to be independent of how the 
compound nucleus is produced. In this spirit one does not need to care for the truly initial 
stages of the reaction. One simply assumes the system to be located initially around the 
ground state minimum Q = Qa of the static energy. However, there still is room for the 
precise definition, as the system may or may not be in (quasi-)equilibrium with respect to 
the fission degree of freedom Q — which for large damping belongs to the slowest modes 
present. Since the corresponding momentum may safely be assumed to relax much faster 
one may also let the system start sharply at Q = Qa with a Maxwell distribution in P. 
In any case, for the current across the barrier it takes some finite time to build up. This 
apparent delay of fission was taken |^ Q as an indication for additional possibilities to 
emit hght particles beyond the measure given by Fn/Fx > F„/Fbw- 

In case that the full process is studied in a time dependent picture both for fission 
as well as for particle emission, as done in the typical Langevin codes [^, such an effect 
is included automatically. Problems arise, however, if one tries to imitate this delay in 
statistical codes which are in use for analyzing experimental results. Such codes are based 
on time independent reaction theory. For this reason it is not obvious how this method 
may be reconciled with the picture of fission delay, the "transient effect". In the present 
note we like to shed some light on this problem by exploiting the concept of a "mean first 
passage time" (MFPT). 

2 Time dependent current across the barrier 

Solutions of the transport equation require appropriate boundary conditions for the coor- 
dinate Q (and P if present). For the solutions discussed above the boundary conditions 
are chosen to make sure that the distribution vanishes at infinity. Calculations of the 
current j{t) across the barrier then typically imply a behavior as exhibited in Fig.|l| for 
different initial conditions. In all cases the asymptotic value of jh{t) is seen to behave like 
Fk exp(— Fk t/^). The differences at short times are due to the following choices: 

^For Kramers' equation proper this involves the coordinate Q and a momentum P. The latter is 
absent in the Smoluchoski equation into which the former turns into for overdamped motion. 



(i) For the dashed curves the system starts out of equihbrium for both Q and P; the 
curve on the left corresponds to the current at the barrier jb(^) = j{Qh,t) and the right 
one to that in the scission region jsc{t) = j{Qsc,t). The equilibrium is defined by the 
oscillator potential which approximates the V{Q) around Qs,- 

(ii) For the fully drawn line the system starts at Qa sharp with a Maxwell distribution 
in P. The obvious delay is essentially due to the relaxation of Q to the quasi-equilibrium 
in the well. This feature is demonstrated on the right by the t-dependence of the width 
in Q (as given by the average potential energy). 

(iii) The dotted curve corresponds to an initial distribution which is also sharp in the 
momentum but centers at the finite value {P)t=o = V MT. 




Figure 1: The current across the barrier for different initial distributions, see text. 

The figure demonstrates clearly a remarkable uncertainty in the very concept of the 
"transient effect", namely that the jh{t) reaches its asymptotic behavior only after some 
finite time rtrans- Even more important appears to be that the whole effect is due to 
the arbitrariness in choosing time zero. If the calculation were repeated at some later 
time to > Ttrans, the same behavior would be seen! This is due to the fact that the whole 
transient effect only comes about because in the initial distribution there is some favorable 
region in phase space from which it is easiest to reach the barrier. This is demonstrated 
in Fig.|^. There, all initial points in phase space are sampled which cross the barrier after 
some time r^. On the right a sufficiently large time was chosen such that most parts of the 
initial distribution have "fissioned". As exhibited on the left, for the much shorter time, 
of the order of rtrans, only a small fraction of points have succeeded to do this, namely 
those which start close to the barrier with a more favorite initial momentum. The vast 
majority of particles are still waiting to complete the same motion but at a later time! 
This feature is very important for several reasons, in particular (i) for an understanding 
of the essentials of the concept of the MFPT 0, ^, P], (ii) that there still is ample time 
for neutrons to be emitted inside the barrier, even for t ^ rtrans- 

The calculations have been performed by simulating the Langevin equations exploiting 
a locally harmonic oscillation, for the following parameters: T = 3 MeV, Ef, = 8 MeV, 
= 1 MeV and rja = I. The potential was constructed from two oscillators, an upright 
one and one upside down, joined at some point between the minimum and the saddle with 
a smooth first derivative. 
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Figure 2: Samples of initial points which overcome the saddle within a given time r^: Left 
part: ~ Tt^ans; right part r, ^ Tk- 

3 The mean first passage time 

For the sake of simplicity we take the example of overdamped motion for which the 
momentum is always in equilibrium such that one only needs to consider the time evolution 
of the fission coordinate Q. The first passage time may then be defined in the following 
way. Suppose that at t = the particle starts at the potential minimum Qa sharp. 
Because of the fluctuating force there will be many trajectories i which will pass a certain 
exit point Qcx once. This process may take the time tj, the first passage time. The mean- 
FPT rjnfpt(<5a Qex) IS defined by the average (tj) over all possibilities. In order to really 
obtain the mean first passage time the i has to be removed from the ensemble once it has 
exited the interval at Qex: the "particle" can be said to be absorbed at Qcx ("absorbing 
barrier"). As the potential V{Q) is assumed to rise to infinity for Q — > — oo, any motion 
to the far left will bounce back: the region Q —oo acts as a "reflecting barrier". 

The MFPT r^fpt can also be calculated from solutions K{Q,t\ Qa,0) of the Smolu- 
chowski equation, adequate for overdamped motion 0, ^, 0. The initial condition for the 
particles to start at Qa is then: limj^o -^(Q, ^IQa, 0) = 6{Q — Qa), which is identical to 
the one used for the fully drawn line of Fig.[l|. It is this case which (for constant friction 
and temperature) allows for an analytic form for the MFPT, 

T"mfpt(Qa Qex) = 7f ^'^ ^Xp 

by reasoning as follows: The probability of finding at time t the particle still inside the 
interval (— oo, Qex) is given by W{Qa, t) = j^""^ dQ K{Q, t \ Qa, 0). Hence, the probability 
for it to leave the region during the time lap from t to t + dt is determined by —dW = 
— ■^W{Qs,,t)dt, such that one has Tmfpt(Qa — ^ Qex) = — / tdW which turns into 

Tmfpt(Qa ^ Qex) = / / K{Q,t\Q^,0)dQdt= dttj{Q,^,t\Q^,0). (3) 

Jo J -oo Jo 

These formulas are associated to the special boundary conditions with respect to the 

coordinate mentioned before, the reflecting barrier at Q — >■ — oo and an absorbing 
barrier at Qex- In particular for the latter reason it is not permitted to use in (^ the 
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Figure 3: MFPT for the cubic potential normalized to its asymptotic value. Solid, dashed, 
dotted-dashed and dotted curves correspond to T j E^, = 0.1, 0.2, 0.5, 1.0. 

currents shown in Fig.|l[ Inserting them blindly would indeed lead to expressions for 
Tmfpt in which the Xtrans appears^. This is in clear distinction to the correct expression 
Actually, the derivation of (0) involves proper solutions of that equation which is 
"adjoint" to the Smoluchowski equation, and which describes motion backward in time. 
In Fig.^ we show the dependence of r^fptiQa. Qcx) on Qcx as given by (|^) calculated 
for a cubic potential. Evidently, the MFPT needed to reach the saddle at Qb is exactly 
half the asymptotic value. The latter may be identified as the mean fission life time 
Tf = T^{pt{Qa — > Qex ^ Qh) ■ For the typical conditions under which Kramers' rate 
formula (|l|) is valid for overdamped motion, this relation of Rk = (t/)^ to the asymptotic 
MFPT can be derived analytically [||. Another remarkable feature seen in Fig.|^ is the 
insensitivity of the MFPT to the exit point for small and large Qex- 

4 Discussion 

It should be evident from the previous discussion that in the very concept of the MFPT 
there is no room for a transient effect. After all, formula (|]) is based on exact solutions of 
the transport equation which satisfy the same initial condition as those used for the plots 
in Fig^. One essential difference is seen in the fact that the evaluation of the MFPT 
takes into account an average over all initial points, as is warranted by the definition of 
the MFPT through the probability distribution —dW. Contrasting this feature, and as 
outlined in sect.|^, the transient effect only represents a minor part of the initial population, 
namely that one which reaches the barrier first. Discarding the rest implies ignoring the 
many particles which are still moving inside the barrier for times typically much longer 
than rjnfpt(Qa Qh)- Neutrons from deformations corresponding to that region are not 
only emitted within Xtrans but within rmfpt(Qa Qh), which turns out to be just half of 
the total fission time tk. Of course, this discussion shows that — besides the additional 
neutrons often associated to the transient effect — also those, supposedly emitted during 
the motion from saddle to scission are not treated correctly by introducing the saddle to 

•^Approximating the form of the current by jb(i) = inoim(l — 6^'"'""'=*/'') exp (— FKi/^) one would 
falsely get r^iptiQh) = fi/i^K + ^ trans 



scission time Tssc of |T31- As one may guess from Fig.^, like the rtrans, the Tssc does not 
appear to be in accord with the MFPT either. These findings suggest that one simply 
estimates the emission rate of neutrons over fission from the ratio Tn/Tm. Anything else 
does not seem to be in accord with the assumption that fission can be described by the 
Kramers or Smoluchowski equations, for the usually assumed form of the potential. This 
does not rule out any effects related to a more complicated dynamics, in particular if the 
initial stage of the whole reaction is to be described with a different transport theory. 
These findings may perhaps imply that some of the existing statistical codes will have to 
be revised. The question of the temperature dependence of nuclear transport does not 
seem to be a closed one yet. A theoretical prediction has been given in 0|. 
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